Structural Determination of Multilayer Graphene via Atomic Moire Interferometry 
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Rotational misalignment of two stacked honeycomb lattices produces a moire pattern that is observable in 
scanning tunneling microscopy as a small modulation of the apparent surface height. This is known from 
experiments on highly-oriented pyrolytic graphite. Here, we observe the combined effect of three-layer moire 
patterns in multilayer graphene grown on SiC (0001). Small-angle rotations between the first and third layer are 
shown to produce a "double-moire" pattern, resulting from the interference of moire patterns from the first three 
layers. These patterns are strongly affected by relative lattice strain between the layers. We model the moire 
patterns as a beat-period of the mismatched reciprocal lattice vectors and show how these patterns can be used to 
determine the relative strain between lattices, in analogy to strain measurement by optical moire interferometry. 



I. INTRODUCTION 

Perfect graphene is a single atomic layer of carbon atoms 
arranged into two interpenetrating triangular sublattices (A 
& B). It has a unique linear band structure stemming from 
the unperturbed n orbitals, which lie above the plane of hy- 
bridized sp 2 bondsPWhen placed on a substrate, the band 
structure may be modified by interaction of the n orbitals with 
the substrate. One way to minimize this interaction is to use 
graphene as its own substrate. In Bernal stacked (AB) bilayer 
graphene, the low-energy band structure near the Fermi en- 
ergy is not linear due to a breaking of the sublattice sym- 
metry; however, if the layers are rotated away from Bernal 
stacking the sublattice symmetry and linear band dispersion 
is preserved over the unit cell of the newly formed moire 
superlatticepl 

In addition to affecting the electronic structure, the moire 
pattern resulting from lattice misalignment or mismatch can 
itself be a useful tool for understanding structural properties. 
Optical moire interferometry has been applied in strain anal- 
ysis for many years.^ Because a scanning tunnneling micro- 
scope (STM) can detect moire patterns in rotated graphene 
layers, 4 analogous methods can be applied to measure local 
strains between individual sheets of graphene at the nanome- 
ter scale. Here we show how STM-based "atomic moire in- 
terferometry" can go beyond surface properties to detect lat- 
tice orientations and strains for depths of several layers. We 
expect this type of analysis to be beneficial to understanding 
electronic and transport properties in graphene multilayers. 

Over the past 20 years, several studies of moire patterns 
occurring on highly-oriented pyrolytic graphite (HOPG) have 
been conducted. 5 10 These typically occur due to a rota- 
tion between two layers near the surface, or from an exfoli- 
ated flake resting on the HOPG surface. 11 Graphene islands 
grown on transition metal surfaces (e.g., Ir(lll), Ru(0001), 
or Pd(l 1 l)jl2lEEl a i so display moire-style superlattices due to 
both rotational misalignment and the inherent size mismatch 
of the lattices. More recently, it was found that multilayer epi- 
taxial graphene (MEG) grown on SiC(OOOl) stacks in such 
a way that almost every pair of graphene sheets is rotated 
with respect to its neighbor(s). Graphitic Bernal stacking 
of the graphene layers appears to be rare in this material, 



although there are clearly preferred layer orientationsPEEH 
X-ray diffraction and low-energy electron diffraction experi- 
ments show broad peaks corresponding to rotations of ±2.2° 
and 30° to the SiC lattice.™ The rotational disorder of 
the graphene samples produces moire patterns in ST M to- 
pographs, and has been under active investigation ! 2 ! 4 ! 16 ! ^ The 
patterns appear as an additional corrugation in STM imaging; 
however, the apparent height variation is thought to be dom- 
inated by a modulation in the local density of statesP122l The 
local stacking changes continuously between regions that re- 
semble AA, AB, BA and slip stacking. 7 Despite regions of lo- 
cally AB stacking, the sublattice symmetry is preserved over 
the full moire unit cell. In each unit cell, AB regions are com- 
plemented by BA regions that restore the inversion symmetry. 
The result is that the electronic properties of an TV-layer film of 
MEG are nearly equivalent to TV independent graphene layers: 
the linear band dispersion near the charge neutrality (Dirac) 
points of the Brillouin zone is preserved! 2 ! 15 ! 21 ! 



H. EXPERIMENTAL METHODS 

Multilayer graphene is grown on SiC (OOOl) by a low- 
vacuum induction furnace technique! 22 ! 23 ! When heated, the 
SiC crystal thermally decomposes leaving behind carbon 
atoms to reform into graphene sheets. The sample studied has 
an average of 10 ± 1 layers as measured by ellipsometry. It 
was initially imaged using atomic force microscopy (AFM) 
and inspected with low energy electron diffraction (LEED) 
and Auger electron spectroscopy (AES) to ensure the sam- 
ple quality. The sample was then placed into a custom built 
ultra-high vacuum chamber (base pressure < 10~ 9 Pa), where 
it was subsequently heated to 1250 °C in order to clean the 
surface after exposure to air. After cleaning, it was studied at 
low temperature (4.3 K) in a cryogenic STM chamber.^ 



III. MULTILAYER MOIRE 

Multilayer graphene often yields more complicated moire 
patterns than HOPG, resulting from the successive rotation of 
layers throughout the depth of the material. Multiple moire 
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FIG. 1: Comparison of two moire patterns in MEG showing multiple layer effects, (a) An image of two similar sized moire patterns interfering, 
resulting in a large unit cell pattern. Layers 1 (Surface Layer) and 2, and 2 and 3 have comparable rotation angles in opposite directions (shown 
in d). (b) zoomed in view of (a), (c) FT of (a), (d) schematic of layer orientation based on the moire pattern in (a), (e) A different superstructure 
that indicates low rotation angles between layers 2 and 3. The spots in the FT relative to the reciprocal lattice vectors were used in this case 
to determine which spots correspond to the rotation between layer 1 and 2. (f) zoomed in view of (e). (g) FT of (e). (h) schematic of layer 
orientation based on the moire pattern in (e). 



patterns are often seen during imaging (Figure [TJ. Unlike 
HOPG, Bemal stacking does not frequently occur; thus, sub- 
surface layers will also have moire alignment, which can be 
observed. Various publications have reported an attenuation 
factor (AF) as socia ted with overlayer coverage of a moire pat- 
tern in HOPGP^A fit of this data yields AF„ = e° 81 ", where 
n is the number of overlayers.^In view of this, a second/third 
layer moire pattern would be visible {AF\ = 2.25) if present 
but the visibility of patterns from deeper layers decays expo- 
nentially. 

The length of a superlattice cell can be evaluated as the 
wavelength of an interference pattern resulting from the mis- 
aligned reciprocal lattice vectors. The difference between a 
reciprocal lattice vector k, and a reciprocal lattice vector ro- 
tated by angle 0, k'. = R{8)kj, is used to determine the length 
(D) and orientation (<p) of the pattern. Since we are interested 
in long wavelength moire patterns (rotation angles, 8 < 5°), 
the only relevant vectors of the interference pattern will be 
Ak = ki - k'. for i = j where i is any of six reciprocal lattice 
vectors. For unstrained lattices, only one pair of vectors k and 
k' needs to be considered because of symmetry. Thus, for two 
rotated lattices of undistorted length, we obtain the equation 
for the size of the moire unit cell length-21251 



D = 



2n 



\k-k'\ 2 sin (0/2)' 
and orientation with respect to the atomic lattice, 



= - - 0/2. 
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Figure [T^ shows an example of two interfering moire patterns 
(the collective interference of three graphene layers), obtained 
by STM measurements. On first glance, one might think that 
the moire pattern constitutes a significant modulation in the 
tip height. In fact, the graphene lattice is extremely smooth, 
such that any small changes in tip height become the foremost 
feature. Height modulations due to moire patterns are typi- 
cally only 0.05 nm peak-to-peak. From Fig. [1J and Fig. [TJ), 
two periodicities can be seen. The large periodicity has a 
length of D « 26.5 nm while the smaller periodicity has a 
length of D w 3nm, still much larger than the atomic lat- 
tice spacing. Because more than a single moire pattern is ob- 
served, we can conclude that this is a combined effect of at 
least three layers (in another area of the sample, 3 coexisting 
moire patterns were found, implying the participation of at 
least 4 graphene layers). From the discrete fast Fourier trans- 
form (FT) in Fig.[lJ, we find that two sets of spots are present 
in reciprocal space. These spots correspond to moire lengths 
of £»! = 2.95 ± 0.05 nm and D 2 = 3.35 ± 0.09 nmP 

Using equations (1) and (2) and a = 0.246 nm as the 
graphene lattice constant, we deduce the rotation angles 
needed to produce a moire of size D\ and D 2 as 9\ = 4.78 ± 
0.07 and 9 2 = 4.21 + 0.11 °. Because A<p = <p x - <p 2 ~ 0°, 
as seen in the Fourier transform, we propose the orientation 
of the top three layers in Fig. [TJl. The first and third layers 
(starting from the surface and counting into the bulk) must be 
nearly aligned in order to obtain the observed A(f> in the FT. 
The measured value is A0 = 0.61 ±0.23 °, slightly higher than 
the expected value of A</> = 0.29 ° from the proposed three 
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layer orientation (Fig. |TJi) and equation (2)P*1 Again, we can 
determine the length of the largest interference pattern from 
Ak using the moire reciprocal lattice vectors from the FT, 



"5 — J ■ \->J 

JD\ + D\- 2D x D 2 cos (A<f>) 

Thus, we calculate a value of D3 = 24.7 ±1.6 nm for the size 
of the largest observed pattern. 26 The uncertainty is due to a 
distortion in the hexagonal shape of the pattern, and is likely 
caused by strain, which will be discussed later. As shown 
in the schematic in Fig. [TJi, the top layer is nearly aligned 
with the third. This creates two patterns of nearly equal size 
and orientation. The angle between the two moire patterns is 
small, producing a large interference pattern. 

In the case of Fig. [T^, we have applied a similar analysis to 
determine the orientation of the top three lattices. The small- 
est k vectors in Fig. [lj were determined to correspond to the 
moire pattern of the first and second layers by comparing the 
moire pattern orientation with respect to the reciprocal lattice 
vectors. We propose that there is a small rotation between the 
second and third layer. The interference of these two moire 
patterns is responsible for the hexagonal set of spots on the 
dotted circle in the FT (Fig. [T^). Consequently, the image 
in Fig. [lj looks qualitatively very different than the one in 
Fig. J. 

As further proof of concept, we image the pattern at mul- 
tiple sample biases for a constant tunneling current. Because 
the imaging is done in a constant current mode, the tip will 
move closer to the surface as the sample bias (Vs) is de- 
creased. This technique is used in imaging Bernal stacked 
graphene bilayers, and shows a transition in the imaging of 
the atomic lattice. 27 At higher biases, every atom is imaged; 
however, as the bias is decreased the tip instead images every 
other atom. 

Figures |2|a-d) show the result of imaging at multiple bi- 
ases on a region within the image shown in Fig. At Vs = 
-0.5 V, a single set of moire maxima is seen, corresponding to 
a moire pattern from the rotation of the first/second layer. As 
the bias decreases, the tip begins to image a second moire pat- 
tern. By Vs = -0.05 V, the tip images only the second moire 
pattern, which apparently results from rotation of layers 2 and 
3. From measured deviations in the expected exponential de- 
cay of tunneling current vs tip-height (/ vs z), we believe that 
by Vs = -0.05 V the tip is in contact with the surface. The dif- 
ference in tip height between V s = -0.5 V and V s = -0.05 V 
is Az = 0.39 nm. The sample biases at Vs = -0.3 V, -0.2 V 
(Fig. |2j?,c) show the transition in imaging the top moire to the 
one underneath. The tip bias, tunneling current and tip shape 
are known to have a large effect on the imaging of a moire 
pattern. 10 This behavior was calculated for a rotated layer on 
top of AB stacked graphite 28 and effects similar to our own 
observations were predicted; however, at present we do not 
understand how the moire patterns of this multilayer system 
depend on tip-sample distance. 




FIG. 2: Multibias imaging on the moire pattern in Figure [Tji. Each 
image is 5 nm x 5 nm, and was taken at the same spatial location and 
tunnel current (100 pA). In (a), only the moire pattern of the top two 
layers is imaged, (b), (c) show the transition as the tip moves closer 
to the surface. In (d), the tip has pushed in far enough to only image 
the lower layer moire pattern. We believe that by V s = -0.05 V, the 
tip is actually in contact with the surface. 



IV. STRAIN 

Relative strain induced between two stacked hexagonal lat- 
tices can have a large effect on observed superlattice wave- 
lengths. This phenomenon is well known in moire interfer- 
ometry as a technique for strain measurement,^ and has been 
observed in large period moire patterns on HOPG. 6 The effect 
appears as a distortion in the real space hexagonal shape of the 
superlattice cell, resulting from unequal distortion of both the 
relative length and angular alignment of the reciprocal lattice 
vectors. Distorted moire patterns are commonly observed in 
HOPG when lattice rotation angles (9) are small. We examine 
the effect of a uniform strain in the top layer relative to an un- 
strained layer below. For simplicity, we demonstrate the effect 
of a relative strain in the arm chair direction on the observed 
superlattice periodicity. 

Strain breaks the symmetry of the reciprocal lattice vec- 
tors, so we now consider the case where the lattice has been 
uniformly strained along the arm chair direction, such that 
\a'.\ = (1 + e,)|«,i, where i = 0, ±. The vector £0 along the di- 
rection of strain does not change angular orientation (cpo = 0), 
so under rotation 

K = -^—RiWo- (4) 
(1 + eg) 

This leads to a change in the wavelength of the interfer- 
ence pattern along the strain direction as given by the law of 




Relative 
\ Uniaxial Stress, 
\ Rotated Lattice 



Equal 
Uniaxial Stress, 
Rotated Lattice 



FIG. 3: Strain induced between two stacked hexagonal (real space) 
lattices in the case where one lattice (black/dashed line) is (a) 
stretched but not rotated and (b) stretched and rotated with respect 
to the undistorted lattice (blue/solid line), (c) displays equal strain 
for both lattices. For rotated lattices under relative strain, the angu- 
lar change in the lattice vectors affect the interference pattern differ- 
ently for each direction. As a result, the moire pattern observed in 
STM will appear stretched or distorted. Significantly higher strain 
is required for similarly sized distortions if the lattices are equally 
strained along the same direction. At the low moire angles (#), the 
lattice vectors (when equally strained) will obtain similarly sized an- 
gular distortions <p t which will nullify each other. 
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FIG. 4: Curves display the length distortion between two vectors 
of the moire cell as a function of the moire rotation angle 8. The 
amount of distortion for a small strain rapidly increases with decreas- 
ing moire rotation angles (i.e. when the rotation angle is comparable 
to the angular distortion, 6 ~ ip). 



In order to approximate the values of e, for a uniform uniaxial 
stress, we make use of Poisson's ratio v and the strain matrix 
e by, e x = e, e v = -vs. 29 The strain matrix is: 



1 + e 
1 -ve 



(8) 



where the axes have been chosen to diagonalize the strain ma- 
trix. For strain in the arm chair direction, we do not need to 
rotate the lattice vectors a, with respect to e. In this case the 
lattice is stressed along ao- We can relate the strain matrix to 
the values e, by applying the original definintion, 

\a\\ = (1 + ei)\ai\ = \eai\. (9) 

Thus, under a small strain approximation (s <K 1) we obtain 

1 



-(1 - 3v)e. 



(10) 



cosines, 



A) 



In 



a(l + eb) 



'*° *°' 2^(l+e )sin 2 (6V2) + e2/4 



Lattice vectors (a ± ) which are not along the direction of strain 
will undergo an additional rotation, <p ± . In the most general 
case of arbitrary changes in length (e,) and arbitrary rotations 
(<pi), this becomes: 



D; 



In 



a(l + e,) 



I*/ - K\ 



2^(1 + e i ) sin 2 t(f5 + Vi )/2] + ef/4 



where i = 0, +, and the distortion in the angular orientation of 
the moire is given by 



arcstn 



(1 + €i)sm(6 + (pi) 



2 A /(l+e ; )sin 2 [(0 + ^)/2] + e 2 /4 / 



fo = 0, tp ± 



: V3i 



1 +V 



4 + 6(1 -3v) 



(11) 



(5) 

We use the values from equations (10) and (11), along with 
v = 0.5 to plot the overall distortion in the wavelength of the 
interference pattern. If we put these values into (6), we can 
compare the length distortion {D + /DJ) of the resulting moire 
pattern. Figure [4] clearly shows the expected behavior. As 
the rotation angle 9 between the lattices increases, the effect 
of the distortion is quickly diminished; however, for a large 
moire pattern the rotation angle is small, leaving the pattern 
(6) easily distorted by the relative strain. The reason is that small 
angular differences become comparable to the rotation angle 
(<p =s 9), and lattice distortions cannot be neglected. Addi- 
tionally, in this regime small changes in the lattice constant a 
contribute greatly to the moire rotation angle. The shape of 
the moire hexagon suffers additional skewing in (pi due to the 
angular distortion of the reciprocal lattice vectors. 

(7) 

Obtaining an exact measure of the strain is difficult because 
it can be hard to determine the direction of the applied stress. 
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FIG. 5: (a) An image of how strain can affect a moire pattern, (b) 
displays the FT of the image in (a). As many as 5 sets of moire 
spots can be identified, requiring interactions from at least 4 layers. 
Strain can be caused by low angle grain boundaries seen in (a) or 
"pleats". The AFM image in (c) shows these structures. The pleats 
are folds in the graphene lattice and tend to run along the SiC step 
or at angles roughly perpendicular to them. The STM image in (d) 
shows atomic imaging on top of a pleat, (e) Schematic of a pleat 
in epitaxial graphene. Mismatches in the thermal contraction after 
growth can induce buckling in the lattice. 

Distortions can also be made from equally stretching both lat- 
tices along one direction as in Fig. [3}:; however the distortions 
of the resulting moire patterns are generally much smaller than 
for relative strains (on the order of several percent strain for 
a comparable length distortion). The reason is that for small 
moire rotation angles, the lattice vectors will obtain similarly 
sized angular distortions for each lattice vector pair, i, thereby 
yielding a relative angular distortion of (pi - <p'. « 0. In this 
case, we can use the original moire equations to approximate 
the length distortion simply by replacing the lattice spacing a 
in equation (1) with its distorted values, a, = (1 + 

Given a low enough angle rotation, the relative strain effects 
are indeed observed. Low angle moire patterns are even more 
likely to occur when multiple patterns can be simultaneously 
observed. Figure[5^ shows a fairly large moire pattern with an 
average unit cell length of D = 10.5 nm + 0.9 nm, correspond- 
ing to a rotation angle 9 = 1.36 + 0.11 °; however, distortion 
from relative strain causes the actual size of the moire pattern 
to range from 9.50 nm to 1 1 .7 nm, depending on which axis is 



measured. This corresponds to a moire length distortion of 
23% (D ma JD mi „ = 11.7/9.50 = 1.23); thus we estimate the 
parameter e w 0.37% for the relative strain between the two 
lattices producing the pattern. 

There are a few mechanisms by which a graphene sheet can 
undergo strain. In Fig. [5^, low angle grain boundaries are vis- 
ible. The moire pattern bends with the boundary on the right. 
In this case, the subsurface rippling of the lattice is likely caus- 
ing a slight relative strain. Another type of strain possible in 
graphene grown on SiC(0001) is a bending of the graphene 
sheet caused by unequal thermal contractions between the SiC 
and graphene sheets after growth!^ In order to relieve stress, 
the lattice buckles (see Fig. |5£). The image in Fig. [5J; was 
taken by atomic force microscopy, and shows these "pleats" 
of the graphene sheet. They are typically 5nm to lOnm in 
height, and are visible in Fig. [3J; as white lines running par- 
allel and transverse to the SiC steps. Despite bending, the 
graphene lattice stays in tact over a pleat, and the graphene 
lattice can extend for several microns before a pleat occurs. 
The STM image in Fig. [5}i was taken on top of a graphene 
pleat. The lattice is continuous, but it is clearly contorted as a 
result of the lattice buckling. The strain reilef looks much like 
the strain relief in nanotubespD Flat areas near these pleats 
could show strain effects. Immediately around pleats, slight 
layer delamination is possible which may not even produce a 
recognizable moire pattern to be observed. 



V. CONCLUSIONS 

Multilayer epitaxial graphene is known to have car- 
rier mea n free times comparable to suspended monolayer 
graphene p 1 * 32 ! 33 ! and measurements to date indicate that its 
electronic structure is effectively that of independent graphene 
monolayers 1151MI However, the physical structure of MEG is 
complex. Using the methods developed here, surface STM 
imaging can be used to obtain layer alignments and local 
strain from at least 3 graphene layers in the multilayer stack. 
These microscopic measurements will be useful for decipher- 
ing the formation of rotational domains and for mapping inho- 
mogeneous strain fields, each of which could affect the quan- 
tum states of this material — especially under the influence of 
electric and magnetic fields. 
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